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Abstract 



We review the problem of the formation of terrestrial planets, with particular emphasis on the 
interaction of dynamical and geochemical models. The lifetime of gas around stars in the process 
of formation is limited to a few million years based on astronomical observations, while isotopic 
dating of meteorites and the Earth-Moon system suggest that perhaps 50-100 million years were 
required for the assembly of the Earth. Therefore, much of the growth of the terrestrial planets in 
our own system is presumed to have taken place under largely gas-free conditions, and the physics 
of terrestrial planet formation is dominated by gravitational interactions and collisions. The ear- 
liest phase of terrestrial-planet formation involve the growth of km-sized or larger planetesimals 
from dust grains, followed by the accumulations of these planetesimals into ~100 lunar- to Mars- 
mass bodies that are initially gravitationally isolated from one-another in a swarm of smaller plan- 
etesimals, but eventually grow to the point of significantly perturbing one-another. The mutual 
perturbations between the embryos, combined with gravitational stirring by Jupiter, lead to orbital 
crossings and collisions that drive the growth to Earth-sized planets on a timescale of 10^ — 10^ 
years. Numerical treatment of this process has focussed on the use of symplectic integrators which 
can rapidy integrate the thousands of gravitationally-interacting bodies necessary to accurately 
model planetary growth. While the general nature of the terrestrial planets-their sizes and orbital 
parameters-seem to be broadly reproduced by the models, there are still some outstanding dynam- 
ical issues. One of these is the presence of an embryo-sized body. Mars, in our system in place of 
the more massive objects that simulations tend to yield. Another is the effect such impacts have 
on the geochemistry of the growing planets; re-equilibration of isotopic ratios of major elements 
during giant impacts (for example) must be considered in comparing the predicted compositions of 
the terrestrial planets with the geochemical data. As the dynamical models become successful in 
reproducing the essential aspects of our own terrestrial planet system, their utility in predicting the 



distribution of terrestrial planet systems around other stars, and interpreting observations of such 
systems, will increase. 
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Dedicated to George Wetherill (1925-2006), pioneer in studies of the formation of the terrestrial 
planets. 
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1 Introduction 



The formation of the terrestrial planets remains one of the enduring problems in planetary science 
and (in view of the expectation of large numbers of extrasolar terrestrial-type planets) astrophysics 
today. The complexity of terrestrial geochemistry, constraints on timescales, the presence of abun- 
dant water on the Earth, and the curious geochemical and dynamical relationships between the 
Earth and the Moon are among the problems that must be addressed by models. Pioneering studies 
by Safronov^ and successors such as Weidenschilling^ established the basic physics of gas-free 
accretion. The effects of gas on accretion were examined somewhat later, most notably by the 
"Kyoto" school of Hayashi and collaborators-. In the 1980's, studies of terrestrial planet forma- 
tion advanced further thanks to George Wetherill"^, his students and postdoctoral collaborators, 
who highlighted the basic problems of obtaining the correct low planetary eccentricities and incli- 
nations, as well as producing a diversity of sizes ranging from Earth through Mars and Mercury. 
Breakthroughs in the subject came through the development of special numerical approaches to 
the problem, as well as theoretical insights that allowed for the right starting boundary conditions. 
Additional geochemical considerations, including formation timescales derived from radioactive 
isotopic ratios, and stable isotopic constraints on source regions, continue to challenge the models 
today. Decades of research have established a rough timeline of events during the formation of the 
Solar System's terrestrial planets. These are summarized in Figure [H which shows the many steps 
which occurred during the formation of the Earth. 

The classical view, developed in the 1960's and 1 970's, is that the planetesimals grow gradually, 
from coUisional coagulation of pebbles and boulders. The growth becomes exponential (runaway) 
when the first massive bodies appear in the disk^-^. However, it is not clear how ordered growth can 
procede beyond 1 meter in size, the so-called meter-size barrier that we explain more extensively 
in section 2. A new view to by-pass the meter-size barrier is that boulders, pebbels and even 
chondrule-size particles can be concentrated in localized structures of a turbulent disk, where they 
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form self-gravitating clumps. The size-equivalent of these clumps can be 10km , 100km or 1,000 
km-, depending on the models and the physics that is accounted for. The growth rate of ~ Moon- 
sized embryos decreases during oligarchic growth because of viscous stirring of planetesimals by 
the embryos and decreased gravitational focusing«i2iii. Late-stage accretion begins when embryo- 
embryo collisions occur-^^, and takes place in the presence of Jupiter and Saturn, which must 
have formed in less than ~ 5 Myr^"*. Late-stage accretion lasted for about 100 Myr in the Solar 
System based on radioisotopic chronometers.—. 

In this review we describe the numerical tools and theoretical concepts used in simulating 
terrestrial planet formation, and the geochemical constraints. We focus on two applications: (1) the 
origin of water on the Earth and (2) the predicted diversity of terrestrial planet systems around other 
stars. We begin by describing the astrophysical and geochemical constraints on timescales. We 
then describe the phases of planetesimal growth and the subsequent oligarchic growth of planetary 
embryos that set the boundary conditions for terrestrial planet formation, following which the 
numerical approach widely used today is outlined. We discuss results from the various groups that 
have conducted simulations, and how well certain constraints from observations are reproduced. 
The relevance of the formation of the Moon by giant impact in understanding terrestrial planet 
formation is considered. We then highlight application of the simulations to the origin of water on 
the Earth, and to simulation of extrasolar planetary systems. We close with a list of outstanding 
issues, and the possible directions for their solution. 

2 Early Phases of Terrestrial Planet Formation 
2.1 Planetesimal Formation 

Planets form from disks created when clumps of interstellar gas and dust, organized in dense 
molecular clouds, collapse to form stars—. The angular momentum content of typical clumps 
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ensures that a portion of the collapsing material ends up in a disk, through which much of the 
mass works its way inward to the growing "protostar" and angular momentum continues to reside 
in the disk-fully consistent with the mass and angular momentum distribution of the Sun and 
planets. Disks undergo evolution from gas-dominated systems to "debris" disks in which only 
solids remain; based on astronomical observations most of the gas is gone within 6 million years 
after the collapse begins (The appearance of the first solids in our solar system is reliably dated 
by meteorites to be 4.568 billion years ago^). 

"Planetesimal" is the term used to connote the fundamental building blocks of the planets 
whose growth is dominated by gravity rather than gas-drag. They are generally defined to be the 
smallest rocky bodies that are decoupled from the gaseous disk. The most commonly-assumed 
planetesimal size is 1 km, corresponding to a mass on the order of 10^^ grams. However, km-sized 
bodies are not completely decoupled from the gas, in that their orbits are significantly altered by 
gas drag via relatively rapid (~ 10'^ — 10"^ yr) damping of their eccentricities and inclinations, and 
much slower (~ 10^ yr) decay of their semimajor axes—. In fact, the actual size distribution of 
bodies during the phase of gravity-dominated growth is determined by the formation mechanism 
of these bodies, which remains uncertain (see below). The planetesimal size is therefore used as a 
parameter in some models of later stages of planetary growth—. 

Modeling planetesimal growth requires a detailed treatment of the structure of the gaseous disk, 
including turbulence, local pressure gradients, magnetic processes, and vortices. Models can be 
constrained by observations of dust populations in disks around young stars, although interpreta- 
tion of observations remains difficult—. There currently exist two qualitatively-different theories 
for planetesimal formation: coUisional growth from smaller bodies?? and local gravitational insta- 
bility of smaller bodies ^i^'^^^^ . 

CoUisional growth of micron-sized grains, especially if they are arranged into fluffy aggregates, 
appears efficient for relatively small particle sizes and impact speeds of ~ 1 m or slowe r ^"^i^^'^^ ; 
see the review by Dominik et al.(2007)^^. However, there is a constant battle between disk tur- 
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bulence, which increases random velocities, and drag-induced settling, which reduces them^^^. 
Growth of particles in such collisions appears effective until they reach roughly 1 cm to 1 m in 
size. At that point, continued growth may be suppressed by collision velocities of > 10 m s^^'^^. 

Meter-sized bodies are the barrier of planetesimal formation. As an object in the gaseous disk 
grows, it becomes less strongly-coupled to the gas such that its orbital velocity transitions between 
the gas velocity, which is slowed by partial pressure support, and the local Keplerian velocity. This 
increases the relative velocity between the object and the local gas such that the object feels a head 
wind which acts to decrease its orbital energy and cause the body to spiral toward the star. Large (> 
10s to 100s of meters) objects have enough inertia that orbital decay occurs slowly, but there exists 
a critical size for which orbital decay is fastest. For the case of rocky bodies in a gaseous disk, 
this critical size is roughly 1 m, and the timescale for infall for meter-sized bodies can be as short 
as 100 years. This is referred to as the 'meter-size "catastrophe" or sometimes "barrier", because 
the infall timescale is far shorter than typical growth timescales^^. CoUisional growth models must 
therefore quickly cross the barrier at meter-sizes if they are to reach planetesimal sizes^^^^^. 

The gravitational instability model for planetesimal formation suggests that a large number 
of small patches of particles could become locally gravitationally unstable and form planetesi- 
mal s^'^'^^ . (The criterion for gravitational instability of Keplerian disks appears already in Safranov 
(1960))^''. This process requires a concentration of meter-sized or smaller particles. If the density 
of solids in a small patch exceeds a critical value, then local gravitational instability can occur, 
leading to top-down formation of planetesimals. A concentration of small particles great by a large 
factor compared with the gas is the key to the process. 

Models for the concentration of small particles often rely on structure within the gaseous com- 
ponent of the disk, generated by turbulence or self-gravity-^^. If the disk is even weakly turbulent, 
a size-dependent concentration of small particles can occur— i^^. Pre-existant chondrule- sized 
particles may have been concentrated at these scales by such a mechanism, thus appearing as the 
basic building blocks of larger structures such as chondritic parent bodies. 
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Self-gravitating clumps of chondrules may end up as 10- to 100-km sized planetesimals; in this 
case particles don't collapse rapidly on the dynamical timescale but slowly contract into planetes- 
imals^. Turbulence can also concentrate larger, meter-sized particles by producing local pressure 
maxima which can act as gathering points for small bodies. As for the meter-size catastrophe, 
boulder-sized objects are the fastest to drift toward pressure maxima-^^13 The concentration in 
these regions can be further increased via a streaming instability between the gas and solids 
and gravitational collapse of the clumps can occur in these dense regions. Johansen et al. (author?) 
[-] showed that planetesimals can form via this process and that the particle clumps (i.e., the rubble- 
pile planetesimals) have a distribution of sizes that ranges up to 1000 km or larger. Figure [2] shows 
the surface density of boulder-sized particles in a disk in which four 1000 km-scale objects have 
formed^. An alternate location for planetesimal formation via gravitational instability are regions 
with an increased local density of solids'*'. Other ways to concentrate solids include drag-induced 
in-spiralling to disk edges (^^, vortices^^, or photo-evaporative depletion of the gas layer—. 

2.2 Oligarchic growth 

Relative velocities in the disk temd to remain low, whether because of damping of eccentricities 
by gas drag coUisional damping, or merely the presence of a few larger bodies that can limit the 
dispersion velocities of the smaller ones. Bodies that are slightly larger than the typical size can 
increase their coUisional cross sections due to gravitational focusing and thereby accelerate their 
growth 

dM _ TlR^T^Vrand ^ j ^ 
dt ~ 2H ^aK^rand), (Ij 

where R represents the body's physical radius, Vrand represents the velocity dispersion of plan- 



4n fact, the idea of the meter-sized catastrophe assumes that the disk has a smooth pressure gradient—. For disks 
with small-scale pressure fluctuations, small particles do not necessarily spiral inward but simply follow the local 
pressure gradient—. 
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etesimals, E is the local surface density of planetesimals, H is the scale height of the planetesimal 
disk, and Fg is the so-called gravitational focussing factor, which depends on Kand- The expres- 
sion for F„ is complicated'*^''*^; for moderate Vrand, it can be approximated by ( 1 + -^^^^ ) , where 
Vesc is the escape speed from the body's surface. 

While random velocities are small, gravitational focusing can increase the growth rates of 
bodies by a factor of hundreds, such that dM/dt ~ M^/^, leading to a phase of rapid "runaway 
growth" ^'^i*°i'*''*"^i'*^i'*^i^° . The length of this phase depends on the timescale for Vrand to increase, 
which depends on a combination of eccentricity growth via interactions with large bodies and 
eccentricity damping. For small (~100 m-sized) planetesimals, gas drag is stronger such that 
runaway growth can be prolonged and embryos may be larger and grow faster-^2i^. 

As large bodies undergo runaway growth, they gravitationally perturb nearby planetesimals. 
The random velocities of planetesimals are therefore increased by the larger bodies in a process 
called "viscous stirring"—. During this time, the random velocities of large bodies are kept small 
via dynamical friction with the swarm of small bodies ' * . As random velocities of planetesimals 
increase, gravitational focusing is reduced, and the growth of large bodies is slowed to the geo- 
metrical accretion limit, such that dM/dt ~ M^/^ (^*'^^). Nonetheless, large bodies continue to 
grow, and jostle each other such that a characteristic spacing of several mutual Hill radii RH,m is 
maintained {RH,m = 0.5[ai + 02] [Mi + M2/3M^]^/'^, where ai and Mi denote the orbital distance 
and mass of object 1, etc.^^). This phase of growth is often referred to as "oligarchic growth", 
as just a few large bodies dominate the dynamics of the system, with reduced growth rates and 
increased interactions between neighboring embryos-^^'^^'^^. 

Figure |4] shows snapshots in time of a simulation of the formation of planetary embryos from 
planetesimals near 1 AU— . Accretion proceeds faster in the inner disk, such that the outer disk 
is still dominated by planetesimals when embryos are fully-formed in the inner disk. Oligarchic 
growth tends to form systems of embryos with roughly comparable masses and separations of 5- 
10 mutual Hill radii-^^^^^. The details of the embryo distribution depend on the total mass and 
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surface density distribution of the disk—. Typical embryo masses in a solar nebula model are a few 
percent of an Earth mass, i.e., roughly lunar to Mars-sized-^^. Figure [3] shows nine distributions 
of embryos with a range in surface densities exponents a and surface densities Si (see Eqn. 1 
of Leinhardt and Richardson ^^). For surface density profiles steeper than r^^, the embryo mass 
decreases with orbital distance. Embryo masses scale roughly linearly with the local disk mass, 
and formation times are much faster for more massive disks. 

The process of embryo formation via runaway and oligarchic growth has very recently come 
into question for three reasons. First, disk turbulence increases the random velocities of planetes- 
imals, often above the critical disruption threshold for km-sized planetesimals. The capacity of 
planetesimals to survive collisions is represented in terms of of Q^, the specific energy required to 
gravitationally disperse half of the object's mass-^^i^. For collisions more energetic than Q}j, col- 
lisions are erosive rather than accretionary making it difficult for embryos to grow. In the presence 
of MRI(magneto-rotational instability) -driven turbulence accretionary growth of large bodies 
appears to require that larger bodies with higher Q}j already exist—. The critical size of these large 
bodies is 300-1000 km. Second, new collision models suggest that planetesimals are weaker than 
previously estimated, such that accretion requires either very slow collisions or pre-seeding of the 
disk with larger objects^"*. Third, statistical models that attempt to reproduce the asteroid belt's ob- 
served size distribution must also resort to seeding the region with large objects of at least 100 km 
in size-^^ . These three lines of evidence all suggest that large, 100-1000 km bodies may have been 
required for the accretionary growth of the much larger embryos. This paradox could be resolved 
if planetesimals form via the turbulent concentration plus gravitational collapse model of Johansen 
et al. (2007) (author?) [^], who inevitably formed 1000 km-scale bodies in MRI-turbulent disks. 
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3 Late-Stage Growth of the Terrestrial Planets 

The planetary embryos formed during the previous oligarchic growth phase begin to perturb one 
another once the local mass in planetesimals and embryos is comparable—. The orbital eccentric- 
ities of embryos become excited, which leads to a phase of close encounters and collisions with 
moderate velocities. Thus begins the final stage of terrestrial planet formation, which ends with the 
formation of a few massive planets. 66,67^68,6970 _ duration of this phase is shortened through the 
presence of Jupiter, which increases the eccentricities of the embryos' orbits and hence the mutual 
collision rates. 

Wetherill (1992)(author?) [— ] was the first to suggest that the formation of planetary embryos 
was not necessarily limited to the terrestrial planet region. He proposed that planetary embryos 
formed also in the asteroid belt. The mutual perturbations among the embryos, combined with the 
perturbations from Jupiter, would have eventually removed all the embryos from the asteroid belt, 
leaving in that region only a fraction of the planetesimal population on dynamically excited orbits. 
For this reason, some of the simulations of Chambers and Wetherill (1998) (author?) [— ] started 
with a population of embryos ranging from ~0.5 to ~4 AU. 

Most recent simulations take advantage of fast symplectic integrators such as Mercury^^ or 
SyMBA^^. These integrators are optimized for planetary studies, and employ algorithms that allow 
for roughly 10 times fewer time steps per orbit as compared with a brute-force N-body integrator, 
for the same accuracy. These integrators also allow for close encounters between bodies, either by 
numerically solving the interaction component of the Hamiltonian (Mercury) or by recursively 
subdividing the time step (SyMBA). When performing integrations with these codes, it is always 
important to choose a time step that is small enough to resolve the orbits of the innermost particles 
with at least ~ 20 time steps per orbit to avoid numerical errors—"^. Collisions are generally 
modeled in a very simplistic fashion, as inelastic mergers occurring anytime two bodies touch. 
Although this assumption appears absurd, it has been shown to have little to no effect on the 
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outcome of accretion simulations However, more complex models show that dynamical friction 
from coUisional debris may play an important role at late stages'-. 

A convenient approximation is often made to reduce the run time needed per simulation, by 
neglecting graviational interactions between planetesimals (see Raymond et al. (2006)for a dis- 
cussion of this issue)^. Assuming that planetesimals do not interact with each other, the run 
time r scales with the number of embryos A'^g and the number of planetesimals, Np, roughly as 
T ~ + 2NeNp. The non-interaction of planetesimals eliminates an additional term. Note 
that T refers to the computing time needed for a given timestep. The total runtime is r integrated 
over all timesteps for all surviving particles. Thus, a key element in the actual runtime of a simula- 
tion is the mean particle lifetime. Configurations with strong external perturbations (e.g., eccentric 
giant planets) tend to run faster because the mean particle lifetime is usually shorter than for con- 
figurations with weak external perturbations. 

Tree codes, which subdivide a group of particles into cells using an opening angle criterion, 
have the advantage over serial codes in that the run time scales with particle number A^ as NlogN 
rather than A^^. They can be run in parallel on several CPUs to further reduce the runtime. Tree 
codes have been used to study planetary dynamics, but to date are only useful in the regime of 
large A^ (A^ > 10'^;'^). The reason for this is that a large amount of computational "overhead" is 
required to build the tree, such that for small A^ more computing time is needed for building the 
tree, and if run in parallel, for communication between processors. The break-even point between 
serial codes and tree codes, for example, is at A^ ~ 1000^°. An advantageous hybrid method for 
large A^ accretion simulations is to integrate particles' orbits with a parallel tree code until A^ drops 
to about 1000, then switch to serial code for the rest of the simulation—. 

A common problem with the current generation simulations is that the final terrestrial planets 
are on orbits that are too eccentric and inclined with respect to the real orbits. The orbital excitation 
is commonly quantified by the normalized angular momentum deficit—: 
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AMD = — ^ (2) 

where aj, Cj, ij, and ruj refer to planet j's semimajor axis, eccentricity, inclination with respect 
to a fiducial plane, and mass. The AMD of the Solar System's terrestrial planets is 0.0018. For 
comparison, the Chambers and Wetherill (1998) Model C simulations, each consisting of at most 
50 bodies extending out to 4 AU and assuming the present orbits of Jupiter and Saturn, yield a 
median AMD of -0.033. The Chambers (2001) simulations 21-24, each consisting of about 150 
bodies and also assuming the present Jupiter and Saturn, have a median AMD of -0.0050. Those 
simulations only extended out to 2 AU, and it is likely that their AMD would be even higher if they 
were extended out to 4 AU (e.g., the Chambers and Wetherill (1998) Model C Simulations, which 
extend out to 4 AU, have a median AMD about 50% larger than in their Model B simulations, 
which only extend to 1.8 AU). 

The missing physics responsible for this mismatch between simulations and constraints is an 
open subject of scientific debate. It has been proposed-^^ that a remnant fraction of the primordial 
nebula would have damped the eccentricities and inclinations of the growing planets. In this case, 
however, the simulations typically form systems of planets that are too numerous and too small. 
The work has been extended by including the effects of secular resonance sweeping as the solar 
nebula dissipates This both forces mergers to reduce the number of final terrestrial planets to 
be comparable to our Solar System, and shortens the growth timescale so that there is sufficient 
nebular gas at the finish to damp the eccentricities to match those of the Solar System terrestrial 
planets. The MHD turbulence of the nebula might also alleviate the problem, enhancing the prob- 
ability that the proto-planets collide with each other and thus leading to systems with a smaller 
number of larger planets—. A problem with both of these scenarios, however, is that since they oc- 
cur on the timescale comparable to the existence of the nebular gas (a few to ~10 Myr). This is not 
consistent with the significantly longer formation timescales inferred from isotopic chronometry 
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of the Earth-Moon system—, discussed more extensively at the end of this section. 

Another possible way to reconcile the simulation results with the constraints is the inclusion 
of dynamical friction. Dynamical friction occurs if embryos and proto-planets evolve among a 
population of small planetesimals with a total mass comparable to the total mass of the embryos. 
A bi-modal distribution of embryos and planetesimals such as this is the likely result of oligarchic 
growth '^^•5'^. Dynamical friction produces the equipartition of the "excitation" energy (e.g., related 
to velocity dispersion, in analogy to the temperature of a gas) between gravitationally interacting 
bodies: the smaller ones obtain higher relative velocities, and the larger ones lower. The relative 
velocity of embryos (hence their eccentricities and inclinations) will therefore be kept low by 
dynamical friction. The simulation of a large number of small planetesimals is, of course, very 
CPU-intensive. Thus simulations typically neglect the effect of the small bodies, or include only a 
limited number of them, which are, therefore, artificially too massive. 

An example evolution of an accretion simulation is shown in Figure[5F^ . This simulation started 
with 1886 sub-embryo sized objects, and is one of the most computationally expensive to date, 
having required 1 .2 x 10"^ CPU hours. The simulation contains a single Jupiter-mass giant planet at 
5.5 AU (not shown), and the evolution is characteristic of simulations with low-eccentricity giant 
planets. Eccentricities are excited in the inner disk by mutual scattering between embryos, and in 
the outer disk via resonant and secular forcing from the giant planet. Dynamical friction acts to 
keep the eccentricities of faster-growing embryos the smallest, and accretion proceeds from the 
inside of the disk outward. Only when embryos reach a critical size can they scatter planetesimals 
and other embryos strongly enough to cause large-scale radial mixing, which is evident in Fig. [5] 
by the change in colors (which represent water contents) of objects. The Earth analog in this 
simulation started to accrete asteroidal water only after ~ 20 Myr of evolution, when it was more 
than half of its final mass. At the end of this simulation, three planets have formed: reasonable 
Venus and Earth analogs at 0.55 and 0.98 AU, and a much-too-massive Mars analog at 1.93 AU. 
Figure [6] shows the growth of the three planets in time. The accretion of the Earth analog occurs 
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on the correct timescale, as it experiences its last giant impact at t 60 Myr. The Venus and Mars 
analogs experience their final giant impacts at 22 and 40 Myr, respectively. 

The comparison between the results in Chambers and Wetherill (1998)(author?) [^^] (no small 
bodies included) and Chambers (2001)(author?) [— ] (accounting for abi-modal mass distribution 
in the initial population) suggested that dynamical friction is indeed important and can drive the 
simulation results in a good direction. Thus O'Brien et al. (2006)(author?) [— ] performed new 
simulations, starting from a system of 25 Mars-mass embryos from 0.5 to 4 AU, embedded in 
a disk of planetesimals with the same total mass and radial extent as the population of embryos, 
modeled with 1 ,000 particles. They performed two sets of four simulations. In one set, called 'EJS' 
for 'Eccentric Jupiter and Saturn,' Jupiter and Saturn are assumed to be initially on their current 
orbits, and in the second set, called 'CJS' for 'Circular Jupiter and Saturn,' they are assumed to be 
on nearly circular orbits with a smaller mutual separation. 

The results of the EJS simulations, with Jupiter and Saturn initially on the current, eccentric 
orbits, can be compared directly to those of Chambers (2001) (author?) [™]. The eccentricities 
and inclinations of the final terrestrial planets measured through the AMD turn out to be five 
times smaller, on average, than in Chambers' runs, and even ~40% lower than in the real Solar 
System. The median time for the last giant impact is ~30 Myr. For comparison, while Chambers 
(2001) (author?) [^°] does not report the time of last giant impact, his Earth and Venus analogues 
take 54 and 62 Myr, respectively, to reach 90% of their final mass. (A recent study by one of 
the authors (SNR) and his colleagues suggest a spread of a factor of a few, and sometimes larger, 
between the last giant impact on Earth analogs in different simulations with the same set of initial 
conditions but different random number initializations). 

The CJS simulations, with Jupiter and Saturn initially on quasi-circular orbits and with smaller 
mutual separations, give a median time for the last giant collision of about 100 Myr. They still 
give terrestrial planets that are a bit too dynamically excited. The eccentricities and inclinations, 
as measured by the AMD, are about 60% larger on average than those of the real terrestrial plan- 
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ets. These somewhat unsatisfactory results with regards to dynamical excitation do not imply 
necessarily that Jupiter and Saturn had to have their current orbits when the process of terrestrial 
planet formation started. The spectacular improvement in the results between the runs in Chambers 
(2001)(author?) [™] and the EJS simulations in O'Brien et al. (2006)(author?) [^] demonstrates 
the dramatic effect of dynamical friction on reducing planetary excitation. With only 1,000 par- 
ticles used to simulate the planetesimal disk, there is no reason to think that the simulations by 
O'Brien et al. (2006)(author?) [— ] give a fully accurate treatment of dynamical friction. Thus, 
it is possible that a future generation of simulations, using more particles of smaller individual 
mass to model the planetesimal disk, and allowing for the regeneration of planetesimals when gi- 
ant impacts occur between embryos^^, would treat dynamical friction more accurately and lead to 
satisfactory results even with Jupiter and Saturn starting on circular orbits. 

Several statistical quantities exist to compare the properties of a system of simulated terrestrial 
planets with the actual inner Solar System—. These include the number and masses of the planets, 
their formation timescales, the AMD of the system, and the radial concentration of the planets 
(the vast majority of the terrestrial planets' mass is concentrated in an annulus between Venus and 
Earth). Reproducing all observed constraints in concert is a major goal of this type of research—. 

With respect to formation timescales, constraints are available from measurement of radioactive 
isotopic systems in rocks on the Earth and Moon. To date these have yielded conflicting results. A 
very detailed analysis^, uses these chronometers, the identity of the tungsten isotopic ratios in the 
Moon and the Earth's mantle, and isotopic dating of the oldest moon rocks. They conclude that 
the last giant impact-that which formed the Earth's Moon-occurred between 50-150 million years 
after the appearance of the first solids in the protoplanetary disk which formed the solar system. 
We assume-but cannot demonstrate- that this giant impact did not occur a significant fraction of 
the Earth-formation time later than the collisions that built the Earth to its present size. With that in 
mind, we argue that any simulations which grow the Earth on a timescale roughly between a few 
tens of millions and 150 million years are consistent with the indications from the geochemical 
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data. 

With respect to the radial distribution of terrestrial planet mass, the simulations described above 
start with a power law column density of solids. In contrast, Chambers and Cassen (2002)(au- 
thor?) [— ] simulated late-stage accretion by generating embryos from the detailed disk model of 
Cassen (2001)(author?) [~] which contains a peak in the surface density at 2 AU (in that model, 
S oc a for a < 2 AU, and S oc a^°-^ for a > 2 AU). They found that simulations with embryos 
generated from a standard MMSN model did a much better job of reproducing the properties of the 
terrestrial planets than the more detailed theoretical disk model. Jin et al. (2008) (author?) [™] 
created a disk model with multiple zones, assuming that the ionization fraction of the gas varied ra- 
dially, thereby affecting the local viscosity and causing pileups and dearths of gas at the boundaries 
between zones. They suggested that non-uniform embryo formation in such a disk could explain 
Mars' small size. Preliminary simulations by one of the authors (S.R.) and colleagues have called 
into question this suggestion. 

3.1 Delivery of Water-Rich Material from the Asteroid Belt 

An outstanding application of the dynamical models is to the problem of the origin of Earth's 
water. The oceanic water content of the Earth is about 0.02% the mass of the Earth, and various 
geochemical estimates put the total amount of water that was present early in the Earth's history at 
5-50 times this number, some or all of which may yet reside in the mantle However, meteoritic 
evidence and theoretical modeling suggest that the protoplanetary disk at 1 AU was too warm at 
the time the gas was present to allow condensation of either water ice or bound water. Therefore, 
there has been a longstanding interest in models that deliver water ice or water-rich silicate bodies 
to the Earth during the latter's formation. Much of the isotopic and dynamical evidence against 
cometary bodies being a primary source, oft quoted in the literature, has been reviewed recently—, 
and a comprehensive treatment of the geochemical evidence is beyond the scope of this review. 



16 



Likewise, alternative models for local delivery of water, for example in the form of adsorbed water 
on nebular silicate grains have been proposed, but will not be described. Of interest here is how 
the dynamical models described above can be used to quantify the delivery of large bodies to the 
Earth from the asteroid belt, where chondritic material (in the form of meteorites) has an average 
D / H ratio close to that of the Earth's oceans. 

The fact that the D / H ratio of Earth's water is chondritic prompted Morbidelli et al. (2000)(au- 
thor?) [— ] to look at dynamical models of terrestrial planet formation to investigate whether a 
sufficiently large amount of mass could be accreted from the asteroid belt. They used simulations 
from Chambers and Wetherill (1998)(author?) [— ], in which planetary embryos beyond 2 AU had 
several times the mass of Mars, and two new simulations with a larger number of individually 
smaller embryos (masses ranging from a lunar mass at ~1 AU to a Mars mass at ~4 AU). They 
found that 1 8 out of the 24 planets formed in the simulations accreted at least one embryo origi- 
nally positioned beyond 2.5 AU. When this happened, at least ~10% of the final planet mass was 
accreted from this source. Assuming that the embryos originally beyond 2.5 AU had a composi- 
tion comparable to that of carbonaceous chondrites (namely with 5 to 10% of mass in water), they 
concluded that these planets would be "wet', i.e. they would start their geochemical evolution with 
a total budget of about 10 ocean masses of water or more. Moreover, Morbidelli et al. (2000)(au- 
thor?) [— ] also studied the evolution of planetesimals under the influence of the embryos. They 
found that planetesimals from the outer asteroid belt also contribute to the delivery of water to the 
forming terrestrial planets, but at a considerably minor level with respect to the embryos. They also 
found that comets from the outer planet region could bring no more than 10% of an ocean mass to 
the Earth, because the collision probability of bodies on cometary orbits with the earth is so low. 
From all these results, they concluded that the accretion of a large amount of water is a stochastic 
process, depending on whether collisions with embryos from the outer asteroid belt occur or not. 
Thus, they envisioned the possibility that in the same planetary system some terrestrial planets are 
wet, and others are water deficient. 
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The findings of Morbidelli et al. (2000)(author?) [— ] have been confirmed in a series of 
subsequent works-22i2L2S. in particular, with simulations starting from a larger number of smaller 
embryos, Raymond et al. (2007)(author?) [— ] concluded that the accretion of a large amount of 
water from the outer asteroid belt is a generic result, and argued that the fact that 1/3 of the planets 
in Morbidelli et al. (2000)(author?) [— ] were dry was an artifact of small number statistics due to 
the limited number of embryos used in those simulations. 

Figure |7] shows the origin of the material incorporated in the final terrestrial planets in the 
O'Brien et al. (2006)(author?) [^-] simulations. The top panel concerns the set of four simulations 
with Jupiter and Saturn initially on circular orbits, and the middle panel to the set with giant planets 
initially on the current orbits. Each line refers to one simulation. Each planet is represented by 
a pie diagram, with size proportional to the planet's diameter and placed at its final semi-major 
axis. The colors in each pie show the contributions of material from the different semimajor-axis 
regions shown on the scale at the bottom of the figure. This represents the feeding zone of each 
planet. The feeding zones are not static, but generally widen and move outward in time^. 

An important difference is immediately apparent between the two sets of simulations. In the 
set with Jupiter and Saturn initially on circular orbits, an important fraction of the mass of all 
terrestrial planets comes from beyond 2.5 AU, and would likely be water-bearing carbonaceous 
material. About 75% of this mass is carried by embryos, the remaining part by planetesimals. 
Thus, the idea that the water comes predominantly from the asteroid belt is supported. However, 
in the set of simulations with giant planets initially on their current, eccentric orbits, none of the 
planets accretes a significant amount material from beyond 2.5 AU. In that case, if the asteroid belt 
is the source of water, it would have to be through objects of ordinary chondritic nature, typical of 
its inner part. We will come back to this idea below. This dramatic difference between the cases 
with eccentric or circular giant planets had already been suggested—, and is explained by several 
authors-^,^, and-^. 

Thus, a crucial question for the origin of the Earth's water is whether it is more reasonable to 
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assume that the giant planets initially had eccentric or circular orbits. The core of a giant planet 
is expected to form on a circular orbit because of strong damping by dynamical friction and tidal 
interactions with the gas diski2i2lii224M. 

Once an isolated giant planet is formed, if the mass is less than about 3 Jupiter masses, its 
interactions with the gas disk should not raise its orbital eccentricit— (but see Goldreich and Sari 
(2003)(author?) ["'^]), but rather damp it out, if it is initially non-zero. In our Solar System, 
however, we don't have an isolated giant planet, but two. The dynamics of the Jupiter-Satum 
pair has been investigate d A typical evolution is that Saturn becomes locked into the 
2:3 resonance with Jupiter. This case is appealing because it may prevent Jupiter from migrating 
rapidly towards the Sun, thus explaining why our Solar System does not have a hot giant planet. 
Even in the case of 2:3 resonance locking, the orbital eccentricity of the giant planets remain small. 
The eccentricity of Jupiter does not exceed 0.007. There are, however, a few cases in which the 
eccentricity of the giant planets can grow—. For instance, if a fast mass accretion is allowed 
onto the planets, the resonance configuration can be broken, and the eccentricity of Jupiter can 
temporarily grow to ~0.1. Also, if the planets are locked into the 3:5 resonance, the eccentricity 
of Jupiter can be raised to 0.035, which is close to its current value. All these cases, however, 
are unstable and temporary, so one has to invoke the disappearance of the disk at the time of the 
excitation, otherwise the planets would find another more stable configuration and the disk would 
damp the eccentricities back to very small values. So, according to our (limited) understanding of 
giant planet formation and gas-disk interactions, a very small orbital eccentricity seems to be more 
plausible, but an eccentric orbit cannot be ruled out with absolute confidence. 

What seems more secure, conversely, is that when the terrestrial planet formation process be- 
gan, the orbits of Jupiter and Saturn had to have a smaller mutual separation than their current 
value. In fact, all simulations agree in showing that the interaction of the giant planets with the 
massive planetesimal disk that would have existed in the early outer Solar System leads to a sig- 
nificant amount of radial migration-i^^JiLJii. in particular, Saturn, Uranus and Neptune migrate 
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outwards, whereas Jupiter migrates inwards. Thus, the orbital separation between Jupiter and 
Saturn grows with time. Recently, a model of the evolution and delayed migration of the giant 
planets has been proposed, and it reproduces fairly well the current architecture of the outer solar 
system-iJ^^. This "Nice" model assumes that Jupiter and Saturn were initially interior to their 
mutual 1:2 mean motion resonance (MMR), and that the orbits of the giant planets at the time 
they cross their 2: 1 MMR were nearly circular. Giant planet migration and the crossing of the 2: 1 
resonance could be delayed for hundreds of Myr, such that the initial configuration would last for 
the entirety of the terrestrial planet formation process—. The second assumption of circular orbits 
at the time of the resonance crossing does not dismiss, in principle, the simulations of terrestrial 
planet formation starting with Jupiter and Saturn on eccentric orbits, because, in these simulations, 
the giant planets eccentricities are damped very fast by the ejection of material from the Solar Sys- 
tem and meet the requirements of Tsiganis et al. (2005)— and Gomes et al. (2005) model after 
a few tens of Myr. One should explain in this case, though, where such eccentricity comes from. 
Conversely, the first assumption of a smaller initial orbital orbital separation of Jupiter and Saturn 
is essential for the success of that model. 

For these reasons, we have performed a new set of four simulations, where Jupiter and Saturn 
are assumed to have initially the current orbital eccentricities, and an orbital separation consistent 
with the Tsiganis et al. (2005) and Gomes et al. (2005) model. The results in terms of final 
eccentricities and inclinations of the terrestrial planets and accretion timescales are intermediate 
between those of the two sets of simulations in O'Brien et al. (2006)(author?) [— ] discussed above 
with AMD values consistent with the Solar System values and a formation timescale consistent 
with the Hf-W age of the Earth-Moon system The origin of the mass accreted by the terrestrial 
planets is presented in the bottom strips of Fig. U\ The planets at or beyond 1 AU, with only one 
exception, receive an important mass contribution from the outer asteroid belt, that is comparable 
to, if not larger, than that from the inner belt. The planets inside 1 AU typically do not receive a 
significant mass contribution from the outer belt, and the contribution from the inner belt is also 

20 



very moderate. 

We believe that we understand, at least at a qualitative level, the differences between the results 
of these new runs and those of the set of O'Brien et al. (2006)*^^ with Jupiter and Saturn on their 
current orbits, in which none of the planets recieved significant outer-belt material. If the orbits 
of the giant planets are closer to each other, the planets precess faster. Thus the positions of the 
secular resonances are shifted outwards. In particular, the powerful uq resonance (occurring when 
a body's perihelion precesses at the same rate as Saturn's), which is currently at the inner border of 
the belt, moves beyond the outer belt. The z/g resonance can drive objects onto orbits with e ~ 1, 
such that they are eliminated by collision with the Sun. It is therefore an obstacle to the transport of 
embryos from the asteroid belt into the terrestrial planet region. In fact, in decreasing their semi- 
major axes from main belt-like values to terrestrial planets-like values, the embryos in the EJS 
simulation have to pass through the resonance. Of course, collisions with the growing terrestrial 
planets are also possible for objects with a Main Belt-like semi-major axis and a large eccentricity, 
but they are less likely. An embryo can be extracted from the resonance by an encounter with 
another embryo, but this is also an event with a moderate probability. So, the flux of material 
from the belt to the terrestrial planet region is enhanced if the z/g resonance is not present. This 
is the case if the eccentricities of the planets are zero as in the CJS simulations (in this case the 
resonance vanishes), or if the planets are closer to each other as in the ECJS simulations (in which 
case the resonance is active, but it is not between the terrestrial planets region and the asteroid 
belt). In the ECJS simulations, the z/g is located around 3.4 AU. We stress that, in order to move 
the uq resonance beyond the asteroid belt, it is not necessary that Jupiter and Saturn are as close 
as postulated in the Tsiganis et al. (2005)— and Gomes et al. (2005)^ model. The initial, less 
extreme, orbital separations used in other models ^'"^^'^ would give a similar result. 

We have recently performed several additional sets of simulations—, including the EEJS ('Extra- 
Eccentric Jupiter and Saturn') set. In four EEJS simulations, Jupiter and Saturn were placed at their 
current semimajor axes but with starting eccentricities of 0. 1 . These systems therefore experienced 
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very strong perturbations from the z/g resonance at 2.1 AU, which acted to remove material from 
the Mars region and also to effectively divide the inner Solar System from the asteroid belt. These 
simulations were the first to produce reasonable Mars analogs, but suffered in terms of water de- 
livery to the Earth. Scattering of embryos and planetesimals during accretion decreased Jupiter 
and Saturn's eccentricities to close to their current values, but the EEJS system does not allow for 
any delayed giant planet migration as may be required by models of the resonant structure of the 
Kuiper belt-^^^^. In fact, it is important to note that the EJS simulations described above are abso- 
lutely inconsistent with the Solar System's architecture because accretion damps the eccentricities 
of Jupiter and Saturn to below their current values, and there is no clear mechanism to increase 
them without affecting their semimajor axes. 

In conclusion, the simulations seem to support, from a dynamical standpoint, the idea of the 
origin of water on Earth from the outer asteroid belt. However, the stochasticity of the terrestrial 
planet accretion process, the limitations of the simulations that we have used, and the uncertainties 
on the initial configuration of the giant planets do not allow us to exclude a priori the possibility 
that the Earth did not receive any contribution from the outer asteroid belt, whereas it accreted 
an important fraction of its mass from the inner belt or its vicinity. For this reason, geochemical 
evidence has been used to try to constrain where the Earth's water came from. For example^ 
have argued that (a) oxygen isotopic differences and (b) siderophile element patterns limit the 
carbonaceous chondritic contribution to 1 % of the mass of the Earth. Constraint (a) can be removed 
or relaxed if the oxygen isotope composition of the Earth and the putative chondritic impactor 
were homogenized in the manner proposed for the Moon-forming impact event (For the Moon- 
forming impactor such a process is deemed essential because the Earth and Moon have identical 
isotopic ratios for both oxygen and tungsten, whereas meteorites vary from these ratios). Constraint 
(b) is a strong one only for relatively small bodies delivering water in a late veneer of material, or 
undifferentiated chondritic embryos mixing fully with the Earth's mantle during the main growth 
phase. If the embryo that delivered the water were differentiated then its core, containing most of 
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the siderophile elements, would not mix with the Earth's mantle. 



4 Extrapolation to Extrasolar Terrestrial Planet Systems 

What counts for terrestrial planet formation? The key parameters are 1) the disk mass and radial 
density distribution, and 2) the giant planet properties (mass, orbit, migration). Here we summarize 
some relevant issues (for a more detailed review): 

• Effect of Disk Properties The accreted planet mass is slightly more than linearly propor- 
tional to the disk mass because the planetary feeding zone widens with disk mass due to 
stronger embryo-embryo scatterin g However, planets that grow to more than a few 
Earth masses during the gaseous disk phase may accrete a thick H/He envelope and be 
"mini-Neptunes" rather than "super Earths" ^^^^'^^ . Whether such objects might be among 
the super-Earth mass planets observed around other stars is an interesting but as yet ill- 
constrained speculation. 

The disk's surface density profile is another key factor. For steeper density profiles, the ter- 
restrial planets form faster and closer to the star, are more massive, more iron-rich and drier 
than planets that form in disks with shallower density profiles—. Disks around other stars 
are observed to have somewhat shallower density slopes^^^^ than the r^^/^ minimum-mass 
solar nebula model of Hayashi (1981)(author?) [— ] and Weidenschilling (1977) (author?) 
[— ]. However, given the preponderance of evidence that giant planets migrate, the validity 
of the minimum-mass solar nebula for either our own solar system or other planetary sys- 
tems is called into question-i^Zii^^. Well-resolved observations of disk surface density profiles 
from facilities like ALMA will help resolve this in the near future. 

• Low-Mass Stars. Low-mass stars are in some sense an ideal place to look for Earth-like 
planets, because an Earth-mass planet in the habitable zone induces a stronger radial velocity 
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signal in the star on a much shorter period than for a Sun-like star— i^^. However, sub- 
mm observations of the outer portions of dusty disks around young stars show a roughly 
linear correlation between disk mass and stellar mass, with a scatter of about 2 orders of 
magnitude in disk mass for a given stellar mass-i^^^^^^ii^. Thus, low-mass stars tend to have 
low-mass disks which should therefore form low-mass giant— and terrestrial planets—. 
However, several low-mass stars are observed to host massive (several Earth-mass), close-in 
planets^^^^. 

• Effect of Giant Planet Properties. Compared with a standard case that includes giant plan- 
ets exterior to the terrestrial planet forming region, the following trends have been noted in 
dynamical simulations: 1) More massive giant planets lead to fewer, more massive terrestrial 
planets 2) More eccentric giant planets lead to fewer, drier, more eccentric terrestrial 
planets 8^i90i97i lOo, 137 _ q[^^i planets have a negative effect on water delivery in virtually all 
cases, overly-perturbing and ejecting much more water-rich asteroidal material than they 
allow to slowly scatter inwards (S. Raymond, unpublished data). 

Hot Jupiter systems represent an interesting situation. In these systems, the giant planet is 
thought to have formed exterior to the terrestrial planet zone, then migrated through that 
zone^. Recent simulations have shown that the giant planet's migration actually induces 
the formation of rocky planets in two ways: 1) interior to the giant planet, material is shep- 
herded by mean motion resonances, leading to the formation of very close-in terrestrial plan- 
g^-g 118,139,140,141,142,143 . 2) extcrior to the giant planet, the orbits of scattered embryos are 
re-circularized by gaseous interactions leading to the formation of a second generation of 
extremely water-rich terrestrial planets at ~ 1 AU-^^^^^. Hence, a key factor is the chronol- 
ogy of migration vs. disk dispersal. If the migration happens when there is still a lot of mass 
in the disk for a good amount of time, then scattered material can be saved and planets can 
formed. 
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5 Conclusion 



Simulation of terrestrial planet formation has become a mature subfield of dynamical astronomy, 
with the potential to provide insight into the origin of our own solar system as well as that of the 
increasing number of multiple planet systems being discovered beyond our solar system. Further 
progress certainly will come from faster computers employing novelties such as, for example, 
many CPUs on a given chip allowing for easy communication between processors and improved 
performance and relevance of parallel codes. But additional insight into the physics and chemistry 
of the problem will be required as well. For example, while the general nature of our terrestrial 
planet system seems to be broadly reproduced by the models, still unexplained is the presence of an 
embryo-sized body, Mars, in place of the more massive objects that the simulations tend to yield. 
Are such outcomes common? We cannot answer this question with the current state of maturity of 
the field. 

Another issue is the effect that collisions between embryos and the growing terrestrial planets 
have on the geochemistry of the latter. The challenge of quantifying in detail the chemical and 
physical processes that occur during giant impacts is a problem outside the scope of the dynamical 
modeling described here, but crucial in trying to relate the geochemistry of the Earth and other 
terrestrial planets to the source material from which they grew. Close collaboration between groups 
that specialize in these two very different types of numerical simulations may permit more detailed 
and confident geochemical predictions in the future. And this, in turn, will increase our confidence 
in the predictions the models described herein can make for the properties of terrestrial planets 
around stars other than our own. 
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Figure 1 : An illustration of the stages of terrestrial planet growth and the relevant timescales 
See 53 for details. 
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Figure 2: Concentration of boulder-sized particles in MRI- turbulent structures in the model of 
Johansen et al. (2007) (author?) [^]. The x and y axes are shown in units of the disk's vertical 
scale height H, and this snapshot is from seven orbital times after a clumping event occurred. 
The greyscale represents the local density of particles, and the solid circles show the location of 
four clumps that are each more massive than Ceres (i.e., they correspond to ~ 1000 km or larger 
"planetesimals" (or small embryos) in the overdense filament. The inset focuses on one clump as 
shown. Original figure provided by Anders Johansen. 



37 



0.2 



0,1 



0.1 







1 1 1 1 1 1 1 1 

- 1 

I a= 1/2 


1 1 1 1 1 1 1 1 
Ej= 10 

a= 1/2 

• "< ; .i 


1 1 1 1 o L 1 1, 1 1- 

Si = . lOD • ° . / . • 
a= 1/2 °'- 

„ O o ' s°°^ 
-J ° ° o =0 

• ■ ' i 

i£Su ° ° '■° 


- 2,= 1 
\ a= 3/2 


Ei= 10 

a= 3/2 ■ 


100 ^ ' ; 

a= 3/2 " 


1 1 1 1 1 1 1 1 1 1 1 1 1 


Ej= 10 ] 

a= 5/2 ■ ■ ] 

lOiiQ 1 1 1 1 1 p 1 '•t 


100 .■ J 

a= 5/2 • • ' . ' 

o 

^ ■ 10 1 



0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 

a[AU] 



Figure 3: Figure 14 from Leinhardt & Richardson (2005) showing their distributions of embryos 
and smaller bodies with a range in surface density exponents a and surface densities Si . All 
panels are at 500,000 years, except for panel 3 in rows 1 and 3 which are at 110,000 and 225,000 
years, respectively. The horizontal bars represent 10 times the Hill radii.This figure was cited 
from Leinhardt & Richardson, Astrophys. J., 625, 427 (2005) with copyright permission from loP 
publishers. 
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Figure 4: Snapshots in orbital eccentricity e vs. semimajor axis a in simulations of the growth of 
planetary embryos The radius of each particle is proportional to the simulation radius but is not 
to scale on the x axis.This figure was cited from Kokubo and Ida, Astrophys. J. 581, 666. (2002) 
with copyright permission from loP publishers. 
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Figure 5: Snapshots in time from a simulation of the late-stage accretion of terrestrial planets, 
starting from 1885 sub-isolation mass objects'*^. The size of each body is proportional to its 
mass^/^, the dark circle represents the relative size of each body's iron core (in the black and 
white version, iron cores are shown only for bodies larger than 0.05 Msarth ), and the color cor- 
responds to its water content (red = dry, blue = 5% water). For a movie of this simulation, go to 
http://casa.colorado.edu/~raymonsn and click on "movies and graphics". 



40 



2.5 F 




0.1 1.0 10.0 100.0 

Time (Myr) 

Figure 6: Growth of the three planets that formed in the simulation from Fig.[5F^, labeled by their 
final orbital distances. The shaded region shows the constraint from isotopic measurements for the 
timing of the Moon-forming impact—. 
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Figure 7: Final terrestrial planets formed in the O'Brien et al. (2006) (author?) [^^] simulations 
(EJS and CJS) as well as the ECJS and EEJS^^ simulations discussed in the text. Pie-diagrams 
show the relative contribution of material from the different semi-major-axis regions, and the di- 
ameter of each symbol is proportional to the diameter of the planet. 
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